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Abstract 

The induction equation induces non trivial correlations between the primordial curvature mode 
and the magnetic mode which is a non linear effect. Assuming a stochastic, gaussian magnetic 
field the resulting power spectra determining the two point cross correlation functions between the 
primordial curvature perturbation and the magnetic energy density contrast as well as the magnetic 
anisotropic stress are calculated approximately. The corresponding numerical solutions are used 
to calculate the angular power spectra determining the temperature anisotropics and polarization 
of the cosmic microwave background, Cg. It is found that the resulting Cg are sub-leading in 
comparison to those generated by the compensated mode for a magnetic field which only redshifts 
with the expansion of the universe. The main focus are scalar modes, however, vector modes will 
also be briefly discussed. 



1 Introduction 

Magnetic fields influence the spectrum of anisotropies of the cosmic microwave background (CMB) 
by contributing to the initial conditions as well as the evolution equations of the perturbations. In 
the following a Friedmann-Robertson- Walker (FRW) background space-time is assumed with a line 
element ds 2 = a 2 (r)(— dr 2 + 5ijdx t dx : >), where a(r) is the corresponding scale factor. Maxwell's 
equations together with Ohm's law in its simplest form, neglecting any convection term, and as- 
suming infinite conductivity leads in a FRW background to PQ-[1], 



(1.1) 



d T (a 2 BJ = V x \v b x (a 2 B 

which yields to [5] 

T (a 2 B) = (a 2 B -v)v b - (v b • v) (a 2 B) - a 2 B ( V • v b ) (1.2) 



where v b is the baryon velocity. Neglecting the righthandside of this equation yields the evolution 
of the magnetic field scaling as a~ 2 . This simplifies the perturbation equations as its contribution 
can be cast in a time-independent way in the mode expansion. This approach has been taken so 
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far in the calculation of the CMB anisotropics [6]-[llj. In most works the magnetic field has been 
assumed to be non helical, except in [11] . As shown first in [12] and investigated in detail in [11] 
helical magnetic fields cause odd parity modes in the cosmic microwave background (CMB) in terms 
of non vanishing cross correlation between the polarization E- and B-modes and the temperature 
and B-mode. 

Equation (jl.2p can be linearized by assuming the magnetic field to be of the form B(x,t) = 
Bq + b(x,r), where Bq is a background field scaling as a~ 2 and b a small perturbation, \b\ <C |5o| 
[2J [13], [3]. Moreover, if in addition it is assumed that on the relevant scales spatial variations of 
Bq can be neglected in comparison to those of b then there is no mode coupling and it is possible 
to identify the three different types of MHD modes. These are the two magnetosonic modes and 
the Alfven mode whose evolution can be studied separately [13] . Whereas the magnetosonic modes 
perturb the density as well as the baryon velocity and thus correspond to scalar perturbations the 
Alfven modes are vector perturbations as the density is left unperturbed. 

In the following the background magnetic field is not assumed to be uniform but rather a 
gaussian stochastic magnetic field, which is assumed to be non helical, for simplicity. The aim is to 
investigate the effect of the evolution of the magnetic field due to the baryon velocity. Assuming 
the presence of a primordial curvature mode generated during inflation it is shown that there is a 
non vanishing correlation between the magnetic field and the inflationary curvature mode. This is 
used to numerically calculate the effect on the CMB. The main focus are scalar modes, however, 
vector modes will also be briefly discussed. 

2 Cross correlation functions between the primordial curvature 
and magnetic mode 

The magnetic field is assumed to be non helical and gaussian so that its two point function is 
assumed to be of the form 



where, for simplicity, a zero correlation time is assumed leading to a Markovian stochastic process 
with a delta function correlation in time. This is in the line of the quasi-normal approximation 
which is used in the problem of closure resulting from an infinite hierarchy of moment equations in 
magnetohydro dynamic turbulence [5]. 

2.1 Scalar modes 

In the longitudinal gauge the fluid velocity coincides with the gauge invariant velocity V. Using 
the expansion in scalar harmonics Q(°\k, x) (e.g. [2]) and using for the magnetic field 



{B*{k,T)Bj{k> \t')) = ^.t^PbCM Ui 




) 



(2.1) 




(2.2) 



k 



then equation (jl.lj) yields to 




(2.3) 
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The baryon velocity is determined by (e.g. [TU], [Z]-[2]) 



V h = (3 C 2 - l)HV h + - 3c 2 ^) + kc 2 s A h + Rt- 1 ^ - V h ) + ^kL, (2.4) 

where, in terms of the magnetic energy density Ab and anisotropic stress ttb (cf. Appendix), 
L = Ab~ §7Tb is due to the Lorentz force J x B and R = f^-- Furthermore, c 2 = g| is the adiabatic 
sound speed and r" 1 is the mean free path of photons between scatterings given in terms of the 
number density of free electrons n e and the Thomson cross section or, t~ 1 = an e aT- Moreover, <& 
and denote the gauge invariant Bardeen potentials and T~L = ^. As can be seen for example in the 
baryon velocity equation, the magnetic field does not enter linearly, but rather quadratically. This 
is also the case for the remaining perturbation equations (e.g., [7]-[ll|). Moreover, the equation 
governing its evolution is an integro-differential equation in the continuum limit. This makes an 
exact treatment of the magnetic field very difficult and therefore here it is assumed that the change 
in the magnetic field due to the effect of the baryon velocity is small so that 

Bi(k,r) = Bf\k) + b l (k, T ), (2.5) 

where \bi\ -C l-B^I- The time- independent solution B- (k) results in constant contributions of the 
magnetic field to the perturbation equations in terms of its energy density contrast and anisotropic 
stress. Moreover, these also enter in the initial conditions [7j- [TT] . The linear perturbation equations 
and initial conditions can be separated in one part proportional to and one proportional to 7Tg. 
The final CMB temperature and polarization angular power spectra are determined by summing 
the two contributions using the auto- and cross correlation functions of and ttb OHO]. Using 
(|2.5p in equation (j2.3j) the back reaction of the fluid dynamics onto the magnetic field is determined 
by 

b t (k, r)=J2 m^±J^L B f\% - q)VM r). (2.6) 

In the standard ACDM model the CMB temperature anisotropies and polarization are due to the 
adiabatic mode which is a primordial curvature perturbation, C(^)> generated during inflation. ( 
is assumed to be a gaussian random variable with two point correlation function (£(/c)*£(A/)) = 

'P({k)Sj:g and the power spectrum is defined by V^(k) = -p-^-s ( -jr) -In the numerical solution 
the best fit values of the 6-parameter ACDM model of WMAP7 will be used [15]. Taking into 
account the contribution from a standard adiabatic mode the total baryon velocity is then given 
by, at zeroth order in the perturbation of the magnetic field, 

V b (k, r) = ((k)V b Hk, r) + A fl(0) (k)V b AB (k, r) + n B{0) (k)V? B (k, r), (2.7) 

where V b , X = (,Ab,ttb, denotes the baryon velocity calculated only considering the contribution 
X. Therefore the first order perturbation of the magnetic field is determined by 

fc&r) = £ m " k ' Bf {k - q) f T dr' [c(g)F 6 C (g, r') + A B(0) {q)V^ B (q, r') 

- Q J Ti 

q 

+W<7)VT0?y)] (2-8) 
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where Tj is some initial time and we assume bj(Tj) = 0. Therefore at first order the magnetic 
field contributions, Ab and ttb, and the primordial curvature perturbation have non vanishing 
correlations. It is interesting to note that in the tight-coupling limit of photons and baryons, long 
before last scattering, equation (|2.6p implies, 



r)^-- A J2 m \' ^ Bf>\k - q)A,(q, r), (2.9) 



using V 7 ~ Vf, and the evolution of the photon density contrast (e.g. [14]) 



A 7 (£,r) = --kV 7 (k,r). (2.10) 
3 



Therefore in the tight-coupling regime the change in the magnetic field is driven by the sum over all 
modes of the photon density perturbation weighted by the amplitude of the initial magnetic field. 
A similar, but simpler relation between the magnetic field perturbation and density perturbation 
was found in [13] for an overall constant mag netic field . 

In order to estimate the effect on the angular power spectra of the CMB the final expressions are 
obtained using As(k,T) and TTB(k,r) calculated by employing (|2.8p . Moreover, relevant definitions 
can be found in the Appendix. This is in contrast to previous work where only the amplitude of the 
magnetic field in /c-space, B^°\k), is taken into account which is constant in time yielding time inde- 
pendent expressions, A B ( )(k) and Tr B (o)(k) [S]-[ll]. As can be seen from (|2.8p the cross correlation 
between the amplitude of the adiabatic, primordial curvature mode £(/c) and the magnetic energy 
density contrast A#(A;,t) as well as the anisotropic stress 7Tb(£;,t) are non zero. Therefore there is 
an additional contribution to the CMB anisotropies. It is assumed that initially, there is no cross 
correlation between the primordial curvature mode and the mag netic mode, (C : (k?)B$ 0) (k)) = 0. 
At some later time r these cross correlations are non vanishing and found to be, at lowest order, 

<C*(*W*,r)> ~ [(C(k')h(q,r)B^\k-q)) + (C(k')Bf\q)b\k - q,r))] (2.11) 

which in the continuum limit together with the expression for the initial magnetic field power 
spectrum 

where ps is the magnetic energy density today calculated using the magnetic field strength smoothed 
over the magnetic diffusion scale k m ~ 200.694 (^) Mpc" 1 with a Gaussian window function 
(cf. [10], [1], [3]) assuming ns > —3, leads to 



r 2 A s p B f^_\ 3+nB fk n 
J W k 3 2T (nB+Z)p^ \k m ) \k p/ 



(C(k')A B (k,r))=S^^^^ (A) / dzz > e • 



l 

< / t dy 



k > 2 - 



e IwJ 2z '^ (i _ 2zy + z 2 ) 2 z [-2(1 + z 2 )y + z(l + 3y 2 )] 
+z nB (-l + 2zy-y 2 )}F vC (x), (2.13) 
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where y = k ■ q/ (kq), z = q/k and 



(x) = f dx'v£(x') (2.14) 



for x = kr and x^ = kr^. Similarly, for the cross correlation between the primordial curvature mode 
and the magnetic anisotropic stress it is found that 

2 o A „_ / u \ n s+ 3 / u \ r°o ( k \ 2 



k \ J2 







dy 



1 - 2zy + z 2 ) 2 (2z 2 + (2 - z 2 )zy - 6z 2 y 2 + 3z 3 y 3 ) 



+z n B (l + zy + y 2_ 3^3)] F ^ (215) 



The two-point correlation function of a random variable % can be written in terms of the dimen- 
sionless spectrum V( xx ) 

2vr 2 

(^Xp) = ^V {xx) (k)6 rk ,. (2.16) 

The dimensionless spectral functions determining the cross correlations between the primordial 
curvature mode and the magnetic mode , "P{^a b )(^i t ) an d 'P(£ 7 r B )(k,T), are shown in figure [T] for 
the best fit values of the 6-parameter ACDM model of WMAP7 [15] at present time, in particular, 
A s = 2.43 x 1(T 9 , n s = 0.963, k p = 0.002 Mpc" 1 , n b = 0.00227/i" 2 , Q A = 0.738, h = 0.714, and 
the reionization optical depth r re = 0.086. Moreover, using a Gaussian window function [10] the 
smoothed magnetic field strength is set to Bq = 10 nG and the spectral index ub = —2.9. These 
values are used throughout in all numerical solutions. The baryon velocity for the primordial 
curvature mode has been calculated using CMBEASY [16J. As can be appreciated from figure 
[T] whereas the primordial curvature perturbation is correlated with the magnetic energy density 
{upper panel) it is anti-correlated with the magnetic anisotropic stress (lower panel) in a large part 
of the (k, r)-plane. Moreover, close to decoupling the acoustic oscillations in the baryon-photon 
fluid manifest themselves in the cross correlations of the curvature mode and the magnetic mode. 
Visually this results in the "ripples" in the color maps (right-hand- side column). The amplitudes 
of the dimensionless spectral functions determining the auto and cross correlation functions of 
constant A B ( ) and ir B (o) for Bq = 10 nG and ub = —2.9 for the same range of k values are of the 
order of 10~ n for (vr^ (0) 7r S (o)) as well as ( A* (0) 7T B (o) ) and of the order of 10~ 12 for (A^ (0) A B ( )) 
|10j . Thus comparing these zeroth order auto and cross correlation functions with the first order 
cross correlation functions in figured] (right-hand- side) around decoupling shows that the latter are 
significantly lower than the former. 

In |17] a non trivial cross-correlation between the primordial curvature mode and the magnetic 
field resulted during inflation by assuming that the magnetic field is generated during inflation and 
conformal invariance is broken by coupling the electromagnetic field to the curvaton field. This is 
different from the case considered here where the non vanishing cross correlations of the magnetic 
field and the primordial curvature perturbation are due to the backreaction of the baryon fluid onto 
the magnetic field. 

Finally, estimating (A* B (k, r)A B (k', r)), {A* B (k, r)-K B (k', r)) and {ir B (k, T)ir B (k', r)) shows that 
the first order corrections to the zeroth order expressions are of the order of [P B (o) (k, k m )] 3 whereas 
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Figure 1: Top panel: A color map of the values of the dimensionless spectral function determining 
the cross correlation between the primordial curvature mode and the magnetic energy density 
contrast and that of the adiabatic mode P((a b )(^ t )' Bottom panel: A color map of the values 
of the dimensionless spectral function determining the cross correlation between the primordial 
curvature mode and the magnetic anisotropic stress and that of the adiabatic mode (k, r). 

In both panels, the figure on the left shows the whole range of r and on the right the time around last 
scattering r/ s = 285 Mpc is shown at a higher resolution. Present time corresponds to tq = 14350 
Mpc. For the primordial curvature mode the bestfit values of WMAP7 are used. 



the zeroth order expressions are of the order of [P B (o)(k,k m )] . Taking into account that the 
magnetic field spectrum is nearly scale-invariant ns = —2.9 these corrections are suppressed, with 
respect to the zeroth order expressions, by a factor ^ ~ 10~ 7 (^) [10J, which for the values at 
hand is of the order of 10 -5 , and will therefore be neglected. This is the same factor which controls 
the corrections to the spectral function determining the two point function of the magnetic field. 



2.2 Vector modes 

Before closing this section we briefly comment on vector modes. In this case the perturbation 
of the fluid 3- velocity is determined by the gauge-invariant amplitude of the matter velocity field 
yf 11 ) (fc, r) [14j . so that the induction equation (jl.2p . implies in terms of the helicity basis, defined 
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in the Appendix (cf. equation (|6.6p ) 

Bi(k, t) = TiYl [( g t)^ ~ ( i t) m km6 ^} B ^ ~ T ) V b ±X \& T )- ( 2 - 17 ) 

Then for Bi(k,r) = + bi(k,r), where |6j(fc,r)| <C |-B- (fc)| the first order solution is given 

by 

h(k,r) = TiJ2 [(e±).g,- - {^)jm5i 3 ] Bf\k-q) vrgj (a) f dr'V^ (q,r f ) (2.18) 

where for the baryon velocity a similar notation to equation (12. 7h was used. The CMB anisotropies 
are determined by the two-point function to first order 

<4 +1) *(fc,T)4 +1 >(£',T) +^~ 1 >(k,T)^ 1 \k',T)) ~ (^*(k)^(k') + ^*(k)^(k')) 

+<^TO ) ^^4 +1) (^^)+^*(^)4" 1) (^r)> + (2-19) 

The scaling with the number of powers of the magnetic field spectrum P B (o)(k,k m ) is the same 
as in the scalar case. Therfore the corrections to the zeroth order expression are suppressed by a 
factor — [11], which for the values at hand is of the order of 10 -5 , and will therefore be neglected. 

3 Resulting angular power spectra 



The line-of-sight integration approach [18] allows to write the brightness function for each compo- 
nent as [19], for the scalar mode, 

@ "J T _lf = jH drS* (k, t)u [k(r - r)} (3.20) 

where S@ is the source function and X denotes and ttb- Thus the angular power spectra 

determining the temperature autocorrelation function are given by 

c TT, m = _i_J™dk £ dTVm{k:T)s C eikiT)s n {kjT) {je [k[To _ r)])2 (3 21) 

where H = Ab,ttb and a similar expression for the autocorrelation functions of the E-mode. The 
cross-correlation between temperature (T) and polarization (E) is determined accordingly by 



V m (k,r)Si(k,T)S%(k,r) 
+V m (k,T)Si(k,T)Si(k,T)} ( je [k(r -r)}) 2 . (3.22) 



The resulting angular power spectra are calculated in a modified version of the numerical code of 
|1U| which is based on CMBEASY [IB]. The results for the bestfit values WMAP7 and magnetic 
field as specified before are shown in figures [20 The presence of a magnetic field before neutrino 
decoupling results in a source term for the evolution equation of the amplitude of the comoving 
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Figure 2: The angular spectrum determining the temperature autocorrelation function . Left: 

C t due to the cross correlation between the primordial curvature mode and the magnetic 

mode. Right: The total angular power spectrum due to the correlated magnetic curvature mode, 
£,TT,{CA B ) _|_ £,tt,{(tt b ) com p ar j son w ith that due to a constant magnetic field, C» ' s<0) s(0) ' > + 

2C ,^(A fl(0)7 r B(0) > + c TT, Km , B(0) ) m 



curvature perturbation. Therefore there is a contribution in addition to the primordial curva- 
ture perturbation generated during inflation. After neutrino decoupling the comoving curvature 
perturbation is constant on superhorizon scales. Moreover, the contributions due to the neutrino 
anisotropic stress and the magnetic anisotropic stress cancel invoking the compensated mode [9]. 
The initial conditions for the numerical solutions are set after neutrino decoupling. Hence the total 
comoving curvature perturbation is given by 

c totai = Cg + c sm (3 23) 

where (5^ = — 3 ^ with /3 = hi^-, and tb and r„ being the time of (instantaneous) 

generation of the magnetic field and the time of neutrino decoupling, respectively [9]. However, as 
the additional contribution to the curvature perturbation is proportional to the amplitude of the 
magnetic anisotropic stress it can be neglected in the calculation of the cross correlation between 
the curvature perturbation and the magnetic field variables. Therefore, the comparison with the 
mode due to a constant magnetic field in /c-space, B^°'(k), in figures [2]|3] is done using the numerical 
solution of the compensated mode. 



4 Conclusions 

The evolution of a primordial magnetic field present since long before decoupling has been calcu- 
lated approximately for scalar and vector modes by solving the induction equation by iteration to 
first order. This effective backreaction of the baryon velocity on the magnetic field causes a non 
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Figure 3: The angular spectrum determining the polarization E-mode autocorrelation function 
Cf E . Left: cf Em) due to the cross correlation between the primordial curvature mode and 
the magnetic mode. Right: The total angular power spectrum due to the correlated magnetic 
curvature mode, (jf E '^^ B ^ + (j EE '^ nB ^ [ n comparison with that due to a constant magnetic field, 
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trivial time evolution of the magnetic field in Fourier space in contrast with the constant time be- 
haviour at zeroth order corresponding in real space to the redshifting of the magnetic field strength 
with the expansion of the universe. Observations of the CMB clearly indicate the presence of 
an adiabatic mode which could be a primordial curvature perturbation generated during inflation 
[13] . Therefore, considering a stochastic, gaussian magnetic field, taking into account the effect 
of the baryon velocity induces non vanishing cross correlations between the primordial curvature 
mode generated during inflation and the scalar magnetic mode. This is what would be expected 
by the nonlinear interaction of the magnetic field with the baryon fluid. Using the zeroth order 
numerical solutions for the source functions in the line of sight integration in addition to the spec- 
tral functions determining the cross correlation functions between the primordial curvature modes 
and the magnetic mode the angular power spectra determining the temperature anisotropics and 
polarization of the CMB are calculated numerically. These are the first order corrections due to 
an evolving magnetic field in Fourier space. The numerical results show that these corrections are 
several orders below the zeroth order results and can therefore be neglected, e.g., in the estimation 
of parameters including a primordial magnetic field using the angular power spectra of the CMB 
[20] . The first order corrections to the auto and cross correlation functions of the magnetic energy 
density and the magnetic anisotropic stress of the scalar magnetic mode have been estimated to be 
much below that of the zeroth order expressions and therefore their imprint on the CMB angular 
power spectra will be even smaller. The same conclusion holds for the first order correction to the 
zeroth order correlation function of the magnetic anisotropic stress of the vector mode. 
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Figure 4: The angular spectrum determining the temperature polarization cross correlation function 
Cj E . Left: cJ E ' m due to the cross correlation between the primordial curvature mode and 
the magnetic mode. Right: The total angular power spectrum due to the correlated magnetic 

curvature mode, C» + C« '^ nB> i n comparison with that due to a constant magnetic field, 

c TE,(A b{0) A b(0) ) + 2C ,Ti?,(A s(0) 7r i3{0) > + c TEM B {0)^ B (Q)) ^ 
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6 Appendix 

The energy density of the magnetic field is written in terms of the gauge invariant magnetic energy 
contrast A B such that 

p B = p^A B (k,T)Q^(k,x), (6.1) 

it 

where Q^°\k, x) denote a set of scalar harmonic functions satisfying (A + k 2 )Q^ = (cf. e.g. 
[14|). The magnetic anisotropic stress is determined by 

Hij) {x,r)= Pl £ Y.^B ) (k,r)Q^\k 1 x), (6.2) 

m=0,±l,±2 £ 

where m = denotes the scalar part and = k~ 2 Q\^ + ^Q^°\ the vector part is determined 
by m = ±1 and Qif 1 ^ = — j% [Qifj^ + Q^u, ) an ^ ^he tensor modes are given by m = ±2 [19]. 
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However, here only the scalar and vector modes will be relevant. Expanding the magnetic field as 
a 2 Bi(x,r) = a 2 ZkBi(k,T)Q^(k,x) yields to 

A B (k, t) = -L V B, t (q, r)B l (k - q, r). (6.3) 
q 

A convenient representation of the scalar and vector harmonic functions in flat space is given by 

Q(°\k,x) = e il * (6.4) 

Q^(k,S)i = ±-^(^±^2)^* (6.5) 
v2 

The (spatial) coordinate system defined by the unit vectors ei, e 2 and e% is chosen such that £3 lies 
in the direction of k, thus £3 = k. Moreover, in the helicity basis 



= -^|(ei ±ie 2 ) (6.6) 



?i ■ eS = —1 and ei • = and 
of the anisotropic stress are found to be (e.g. 



so that eL • e5 = —1 and eL • e- = and e ? • = 0. With this choice the scalar and vector parts 



2p 7 o 



£ - 9, r)<f B,(g, r) (tf - g*) - £ - g, r^g*, r) 



k 2 ' 



p 7 o y 



e?Y B % {k - q, T)Bj(q, t)P + (ef Y B,(q, r)Bi(k - q, r)k l 



(6.7) 
(6.8) 



For simplicity, 71^ is denoted as ttq. 
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